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Null boundary controllability of a 1-dimensional heat equation with 

an internal point mass 

Scott W Hansen * and Jose de Jesus Martinez * 


Abstract 

We consider a linear hybrid system composed by two rods of equal length connected by a point mass. 
We show that the system is null controllable with Dirichlet and Neumann controls. The results are based 
on a careful spectral spectral analysis together with the moment method. 


1 Introduction 


In this article we prove the boundary null controllability of the temperature of a linear hybrid system 
consisting of two wires or rods connected by a point mass. More precisely, we consider the following system: 


ii — u" = 0, 
v-v" = 0, 

< i = v'{t, 0) — u'{t, 0), 
u{t, 0) = v{t, 0) = z{t), 
= 0 , 


t > 0, X € uji 
t > 0, X € UJ2 
t > 0 
t > 0 


(- 1 , 0 ) 

( 0 , 1 ) 


with either Dirichlet control 


( 1 ) 


= t>0 (2) 

or Neumann control 


H(t,l) = /(t), t>0. (3) 

In the above and throughout this article, ' denotes spatial derivatives and ’ denotes temporal derivatives. 
In addition, u = u{t,x) and v = v{t,x) denote the temperature on uii and UJ 2 , and z = z(t) denotes the 
temperature of the point mass. The initial conditions at time t = 0 are given by 

{ u°(a:) = ii(0, x), X G uji 
V^{x) = v{0,x), XGW 2 
zO = z(0), 

where the triple v^, z°} will be given in an appropriately defined function space. 

System HI) with the homogenous boundary condition 

u(<,l) =0, t>0 (4) 

can be viewed as the limit of the following “epsilon” system with unit density on (—1,1) \ (—e, e) and with 
density I/2e on (—e, e): 

{ Me — Me” =0, t > 0, X G (—1, —c) 

Me — Me" = 0, t > 0, X G (e, 1) (5) 

^Ze - Ze" =0, t > 0, X G (-e, e) 
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where ite, and satisfy the conditions 

-e) = z^{t, -e), z^{t, e) = v^[t, e), 

Ue'{t, -e) = zj{t, -e), zj{t, e) = Ve{t, e), (6) 

- 1 ) = v^{t, 1 ) = 0 , 

for t > 0. In fact, in [8] the authors have shown that under appropriate assumptions of the initial data, 
solutions of dSI) with ([S|) converge weakly to solutions of ([T|) and (|3]) . 

The hybrid system © is a variant of previously studied hybrid models for systems of strings and beams 
with interior point masses. Hansen and Zuazua used the method of characteristics in m to prove the 
boundary null controllability of an analogous string system with an interior point mass. In |13] Littman and 
Taylor use transform methods to prove boundary feedback stabilization of the string mass system. In [1] 
and [2], Castro and Zuazua used method of non-harmonic Fourier series to prove boundary controllability 
of systems of either Rayleigh or Euler-Bernoulli beams with interior point masses. We refer to m, m. 0, 
[18] . [7] and [5] for related results on control and stabilization of systems of beams with end masses. 

Our main results are the following. 

Proposition 1.0.1. System ([1]) with either Dirichlet control or Neumann control @ is null controllable 
in any time T > 0. More precisely, given T > 0 there is a control f G L^{0,T) such that given initial data 
G L^(uji) X L^(uj 2 ) X R we have that {u(T,x),v(T,x), z{T)} = {0,0,0}. 

The solutions in Pronosition ll.O.ll are defined by transposition in the spaces C(0,T- X_i/ 2 ) for the case 
of Dirichlet control and C{0,T;'H) for the case of Neumann control; see Section [3] 

Our general approach is to reduce the control problem to a moment problem. We consider the case of 
Dirichlet control and Neumann control separately in Section [3] 

2 Preliminaries 

We begin with a discussion of well-posedness of the system © with either homogeneous Dirichlet boundary 
condition ([4]) or Neumann boundary condition 

u'(t,l)=0, f>0. (7) 

Given u, v and z defined on uii, UJ 2 and R. respectively, define y = {u,v,z,Y where * denotes transposition. 
Let 


n = L^iuJi) X L^{uj2) X R 


equipped with the norm 

\\y\\l = Uu,v,z)\\l = \\u\\l^ + \\v\\l^ + \zY 
where || • ||(^i is the usual norm in L‘^{uji) for i = 1, 2. In the Dirichlet case (|1]), let 

-ffuji = {uG I u(-l) = 0} 

'&UJ2 = (v e I w(l) = 0} (8) 

d = {{u,v) G til X -§2 I u(0) = u(0)} 

equipped with the norms 

NIL. = * = 1,2 

||(u,z;)||^=|M|^^^+||u||i^. 

One can see that d is algebraically and topologically equivalent to although it will be more convenient 

to think of d as a subspace of di x 1)2 ■ The space 

W = {(u, w, z) G d X R I u(0) = u(0) = z} 
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is a closed subspace of i? x IR with norm ||(u,w,z)||^ = ||(u,?^)||^- In the Neumann case (O, replace the 
definition of 'du )2 in ® by 

(9) 

and otherwise the space W is defined the same way. In either case, it is easy to show (see [S]) that the space 
W is densely and continuously embedded in the space H. Define the operator A : D{A) C 'H —'H by 

[ (f 0 0 \ 

A = [ 0 cP 0 (10) 

y—(5od Sod oj 

where d denotes the (distributional) derivative operator, So denotes the Dirac delta function with mass at 
X = 0, and the domain D{A) of A is given in the Dirichlet case (|3]) by 

D{A) = {y€W : u€ v G ^^(wa)}. (11) 

and in the Neumann case © by 

D{A) = {y€W : u€ v G iuj 2 ), v'(1) = 0}. (12) 

When D{A) is endowed with the graph-norm topology 

\\y\\liA) = Mn + \\Ayrn 

it becomes a Hilbert space with continuous embedding in 'H. We can therefore write the homogeneous 
point-mass systems o, o and o, © as 

y(t)=Ay{t), y{0)=y°, t>0 (13) 


where y^ = 

Proposition 2.0.2. The unbounded operator A given by pop in domain D{A) as in dill) is a bijective, 
self-adjoint and dissipative operator with a compact inverse. Furthermore, A is the infinitesimal generator 
of a strongly continuous, compact and analytic semigroup (Tt)i>o. 

Refer to [S] for a detailed proof of the above proposition for the Dirichlet case o, © . As a consequence 
of Proposition 12.0.21 given initial data G "H there exists a unique solution 

y G C'([0,oo);'H) 

to the Cauchy problem (fT^ . If in addition, y° G D{A) then y G (^([O, oo), D(A)). 

In the next subsection it is shown that A has only negative eigenvalues, hence —A is positive, self- 
adjoint it provides an isomorphism: D{A) —>■ Ti.. Moreover, fractional powers of —A are well-defined. Let 
Xi = D{A) and for a G [0,1], define = D((—A)“) and = X'^, the dual space relative to the pivot 
space j-L = Xo of Correspondingly, the semigroup T remains an analytic semigroup on the invariant 
subspaces Xa, 0 < a < 1, and extends continuously to an analytic semigroup on spaces Xa, — 1 < a < 0; 
see e.g., [16] for full explanation. The norm on Xa is given by ||?/||^ = {{—A)°‘y, {—A)°‘y)-H- In particular, 
Xi/2 is the completion of Xi with respect to the norm 

\\y\\l /2 = (--42/,y)o- 

Integration by parts gives 

\\y\\l /2 = {y^y)w- 

Thus, Ai /2 is topologically equivalent to Llo(D) in the Dirichlet case © and {/ G H^{n) : /(—I) = 0} in 
the Neumann case ©. 
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2.1 Spectral analysis for Dirichlet case ([T]), (jl]) 

By ProDOsition l2.0.2l the spectrum a {A) of A is contained in the negative real axis and consists of eigenvalues 
{A„} tending to negative infinity with corresponding eigenvectors {(/JnlneN forming an orthogonal system 
for n. 

Proposition 2.1.1. The eigenvalues {A„}„gN of A in the Dirichlet case are distinct and given by 

A 2 fe = -(fc7r)^, X 2 k-i =-tLl forkeN 
where fik is the k-th positive root of the characteristic eguation 

fj, = 2 cot fj,. (14) 


The corresponding eigenvectors are given by 

( sin(fc7ra;)\ / sm{{l + x) p,k)\ 

sm{kTTx) , ip 2 k-i(x) = sin((l - x)fik) 

0 / \ sinifik) ) 

and tpn € D(A) for all n € N. 

Proof. Look for nontrivial functions (/?„ = {Un,Vn, ZA)* € D{A) such that Ag^n = ^nTn- We use an even 
index in the case that = 0 and an odd index when Z„ ^ 0. The eigensystem corresponding to Z 2 k = 0 
reduces to the problem of finding {U 2 k, V 2 fc) such that 

U2ki^) = ^2kU2k{x), X GUJi 
^2ki^) = ^2kV2kix), X €UJ2 

< U',i0) = P4(0) 

U2k{0) = V2ki0)=0 

y2k{-i) = v2k{i) = o. 

It is easy to check that (p 2 k satisfies the above with X 2 k = —{kir)'^. 

Now consider the case that ^ 2 fc-i 7^ 0. The eigenvalue problem reduces to the problem of finding 
functions {U 2 k-i, l 2 fe-i), and real value ^ 2^-1 such that 

U2k-i{x) = -yU2k-iix), XGWi 

yi'k-iy =-l4^-2k-i{x), xeuj2 

< V^k-iiO) - U'2k-M = -yZ2k-i (15) 

t/2fe-l(0) = V2fe-l(0) = ^2fc-l 
tA2fe-l(—1) = V2fe-l(l) = 0. 

From the boundary condition t/ 2 fe-i(—1) = 14fe-i(l) = 0, we have that the solution is of the form 

U 2 k-iix) = sin ({x + l)tik) 

V 2 k-iix) = Csin {(x - l)pk) 

for some constant C to be determined. The continuity condition C/2/c-i(0) = V2fc-i(0) = ^ 2 /c-i gives 

Z 2 k-i = sin(^fc) = -C'sin(^fc). 

Since Z 2 k-i is nonzero we have that pk is not a multiple of tt. Furthermore, we find that C = —1. Then 
from the third equation in (1151) we see that 


2 cot {pk) = P-k- 


(16) 
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Hence the solution to the eigensystem (ITS|) is 


/U2k-i{x)\ /sm{{l+x)fik)\ 

V 2 k-iix) = sin((l - . 

\ Z2k-i ) \ sin(^fc) ) 

Finally, note that since the function F{^) = 2cot/x — ^ decreases monotonically from +oo to —oo over 
the interval ((fc — l)7r, kir) for all fc G N, there is exactly one root of F in each interval {{k — l)7r, kn) for all 
fc G N. Hence the eigenvalues 


{-(fc7r)^}fcgN U 


are distinct. □ 

Proposition 2.1.2. The sequence {^k} in the Dirichlet case 0 satisfies the asymptotic estimate 

Mfe = (fc-l)7r+^+ C> . (17) 

Consequently, consecutive eigenvalues of A in (fOl) satisfy the gap condition: 

|A„+i-A„|>4 + o(^i^. (18) 

Moreover, the eigenfunctions are asymptotically normalized in the sense that 


lim ||i^„|| = 1. 


Proof. From the end of the previous proof, /r/c = {k— l)7r + efe, where 0 < < tt. The characteristic equation 

m can be rewritten as 


(fc - l)7r + Ck 


= cot Ck 


and thus by monotonicity, 

{k — l)7r/2 < cote/c < fc7r/2. 
Taking inverse cotangent of each term gives 

2 2 
arctan -— < Cfc < arctan — --. 

KTT (fe — l)7r 

Hence by Taylor’s formula we obtain (HZl. 

The estimate (ITS)) can be obtained from 


|A 2 fe+i - A 2 fe| = (fJ.k +1 + k7r)(fj,k+i - kn) 


= ^2A:7r + O 
= 4 + 0 


kiT 


O 


ir 


Finally, it is easy to check that ||+ 2 fc|| = 1 for all fc G N and using estimate (ITTl) that ||+ 2 fc-i|P = 1 + 
C>(fc-2). □ 
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2.2 Spectral analysis for Neumann case ([T]), ([7]) 

As in Subsection 12.11 the eigenvalues of A (denoted A„) form a discrete sequence of negative numbers tending 
to negative infinity with corresponding eigenvectors tpn which form an orthogonal system for 'H. 

Proposition 2.2.1. The eigenvalues {A„}„gN of A in the Neumann case & are given by A„ = — where 
{MnjnGN orc the roots of the characteristic equation 

/r = 2cot2^. (19) 

The corresponding eigenvectors are given by 

( sui{p, 2 k-i{x + 1 )) 
tan/r 2 fc-i cos(^ 2 fc-i(x 
sin ^ 2 fc-i 

( cot/i2fc sin{p.2k{x + 1)) 
cos(^ 2 fc(a; - 1)) 
cos pi2k 

and ifn S D{A) for all n S N. 

Proof. The eigenvalue problem A(pn = ^nTn with = {Un, 14,, G D{A) is the following system: 

U'^{x) = XnUn{x), X&UJi 

V"{x) = A„t4(a:), X &U )2 

< 14(0) - C/;(0) = A„Z„ (21) 

4„(0) = 14(0) = 

4/„(-l) = l7(l)=0. 

First note that the possibility of = 0 leads to the trivial solution. Hence 0 for all n G N. Then 

from the first two equations and the boundary conditions we find that 

Un{x) = SmiUnix + 1 )) 

14 (a:) = C'cos(^„(a: - 1)) 

for some nonzero constant C to be determined. The continuity condition t/„(0) = 14(0) gives 

sin p,n = C cos fin 

and since is nonzero for all n G N we have that C = tan/i„. Then from the third equation in (1211) we see 
that 



pin = - tan/i„ + cot pin 


which is equivalent to the characteristic equation (HU). Hence the corresponding sequence of eigenvectors is 


( sin(/i„(l +x)) 
ianpri cos(/x„(x - 1)) 
sin 


which agrees with (1201) after multiplying by normalizing factors \/2 for n = 2fc — 1 and cot pin for 
n = 2k. □ 

Following the ideas of Proposition 12.1.21 one can prove the following result. 
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Proposition 2.2.2. The sequence {^k} in the Neumann case & satisfies the asymptotic estimate 



(22) 


Consequently, consecutive eigenvalues of A in (USD satisfy the gap condition: 



(23) 


Moreover, the eigenfunctions are asymptotically normalized in the sense that 

lim ||v?„|| = 1. 

n—foo 

3 Proof of Controllability results 

We begin with the case of Neumann control: o, ®. 

3.1 Neumann control 

The dual observation problem to ([T]), ([3]) is 


—u — u" = 0, t > 0 

—V — v" = 0, t > 0 

—z = v'{t,0) — u'{t,0), t>0 

u{t, 0) = v{t, 0) = z{t), t > 0 
u(t,—l) = v'(t, 1) = 0, t > 0 


t > 0, X G uji 
t > 0, X G UJ2 


with terminal data at t = T given by 



lA'ix) = u{T,x), X G oJi 
hF{x) = v{T, x), X G UJ 2 


By letting y = {u, v, z)*, the above problem can be written as 

-y = Fy, y(T) =y^ gU, t > 0. 
Then y G C'([0,T],'H) and is given by 

y(t) = T(T - t)y^; 0<t<T. 


(24) 


(25) 


Let y he a smooth solution of the control problem with smooth / G Lf{{),T). Formal integration by parts 
then shows 



0 = 


/ (u — u")u dxdt + {v — v")v dxdt 



Equivalently, 



(26) 
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Since the functional i{y) := {>(1) is continuous on X 1/2 = W it follows from Propositions 5.1.3 and 10.2.1 in 
m that for solutions of (l24l) there exists C > 0 for which 


\\xit, l)||L 2 (o.r) < CWy’^Wu 

y f Gn. 

(27) 

Hence, equation (1551) uniquely defines y{T) as an element of TL. 

Applying this definition for s G [0,r] 

we see 

yGC{[0,T],n) 


(28) 

and there exists C > 0 for which 



\\y\\L<=^io,T-,'H) < ciWy^Wn + 1/ Il2(o,t))- 

(29) 


As before we have the following lemma. 

Lemma 3.1.1. The control problem ([T]), ([3]) is null controllable in time T > 0 if and only if, for any y^ G TL 
there is f G L'^(0,T) such that 


{y°,TTy'^)n = - [ f{t)v{t,l)dt (30) 

Jo 

holds for all ip" G TL, where y is the solution to the observation problem (1241) . 

Proof. First assume that (1501) holds for all y'^ G TL. Then by (1001) . y{T) = 0. Conversely, if / is a control for 
which y{T) = 0, then (l30l) follows from equation (l26l) . □ 

We are now ready to reduce the control problem o, © to a moment problem. Any initial data 
y^ = v^, z^y in TL for the control problem can be expressed in terms of the eigenfunctions as 

y° = Y1 ynPn (31) 

neN 

where the Fourier coefficients {?/°}„gN belong to fi. Let ijn = {un, Vn, ZnY be the eigensolution of (1551) given 
by 

yn{t,x) = e^’^^'^~*'>ipnix). (32) 


In particular, note that 


Unit, 1 ) = 


y/2e^"-^'^ *)tan^„, n odd 


y^gAn(T t) ^ even. 

Applying these solutions to equation (l30l) we obtain the following moment problem: 

= [ /(T-r)e^"^dT, uGN 

Jo 


^XnT 


(33) 


where 


bn — 


— tan/i„ n is odd 

— 1, n is even 


and by Proposition 12.2.21 a„ = G P'. In particular note that for n = 2fc — I 


kiT 


+ o{k-^) = — + oik-p 


ki: 


(34) 


tan/i 2 fe-i = tan 




and furthermore since tan ^ 2 k-i ^ 0 for all k G N, there exists e > 0 such that 

K\ > -, Vn e N. 

n 

From our estimates of /j,„, A„, and a„, it is easy to show that there are constants K,S > 0 such that 




< Ke 


-Sn-‘ 


n e N. 


(35) 


From equations (1^ and (1^ we see that the series ^ 1/A„ converges, and that there exists a constant p > 0 
such that |Afc+i — Afc| > p for all k G N. This implies the existence of a biorthogonal sequence {6ij(r)}jgN 
(see 0 ) such that 


[ dr = = i ^ 

Jo ( 0 , 

By the method of Russell and Fattorini in [15] we have that there are Mi, M 2 > 0 such that 

IIdill < 


(36) 


It is easy to see that the above implies the convergence of 

/(T-r)=^^e"^^d,(r) 

ieN -1 

which provides a solution to the moment problem (gZI). The proof of Proposition 11.0.II for the case of 
Neuman control ([Sj), is a direct consequence of Lemma [3.1.11 and the existence of the biorthogonal sequence 
{di(r)}igN- 


3.2 Dirichlet Control 

The dual observation problem to (|T]), (0) is 


< 


—u — 

u" 

= 0, 


t > 

0, 

X G uji 

— V — 

v" 

= 0, 


t > 

0, 

X GOJ 2 

—z = 


't,0) - u' 

(i,0), 

t > 

0 


II 

0 ^ 

II 

^ z{t), 

t > 

0 


u{t,- 

■1) 

= w(C 1) 

= 0, 

t > 

0 



with terminal data at t = T given by 


{ v?"{x) = u(T,x), X G uJi 
v'^ {x) = v{T,x), XGW2 

5^ = ziTl 


(37) 


(38) 


and observation Y(t) = By letting y = {u,v,zy, the above problem can be written as a Cauchy 

problem as 

-y = Aly, y(r) = t> 0. (39) 

If y^ € Xi /2 = W then y S (^([O, T], X 1 / 2 ) is given by 

y(t) =T(T-t)y^; 0 < t < T. (40) 
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Let y be a smooth solution of the control problem with smooth / S L^(0,r) and let y be solution of the 
dual problem (l3^ . Integration by parts as earlier results in the identity 

{y{T),y^) = {y°,TTy^)n- [ dt (41) 

Jo 

where (•,•) denotes the duality pairing in X_i /2 x X^j^- 
In the case of the heat equation 

{ q = q" 0 < a; < 1, t > 0 

q{t, 0) = q{t, 1) = 0 t > 0 
g(0,a:) =y° e iLo^(0,l) 0 < x < 1 

it is well known (e.g. a) that for each T > 0 there exists C > 0 for which 

i)I1l2(o,t) < c'lk°llffi(o.i)- 

One can verify that the same estimate holds for solutions of (IMl) in the sense that there exists C > 0 for 
which 

1)IU2(o,t) < C'||y"^||i/2 V G Xi/2. (42) 

Since the semigroup T is strongly continuous on X 1/2 = Hq{VI) it follows that the identity (HD) defines the 
value y(s) for all s G [0,T] as an element of X_i /2 for which there exists C > 0 such that 

l|y||L”(0,T;Jf_i/2) < C'(l|y°l|w + II/IIl2(o,t)) 

and moreover 


yGC'([ 0 ,r],X_i/ 2 ). (43) 

The above estimate (l43l) is sometimes referred to as admissibility of the boundary control operator corre¬ 
sponding to Dirichlet control, and can also be derived in the framework of “well posed boundary control 
systems”: see of [TOl Prop. 10.7.1]. 

Analogous to Lemma 13.1.11 the following lemma characterizes the problem of null controllability of ([1]) , 
(I2|) in terms of the solution y of the observation problem (1551) . 

Lemma 3.2.1. The control problem ([T|), (I5|) is null controllable in time T > 0 if and only if, for any y^ G Tl 
there is f G L‘^{0,T) such that 

{y°,'^Ty'^)H= [ f{t)vit,l)dt (44) 

holds for all € 7^, where y = (w, {?, zY is a solution of dMl). 

We are now ready to reduce the control problem (HD, (HD to a moment problem. Any initial data 
yO = (yO^ ^oy control problem can be expressed in terms of the eigenfunctions as 

2/° = XI 

new 

where the Fourier coefficients {y°}rt6N belong to P. Let y„ = (u„, u„, 5„)* be the eigensolution of (1551) given 
by 

yn(t,x) = e^''^'^~*'^ipnix). (46) 

In particular, note that 

, , J1)^, n = 2k 

= n = 2k-l. 
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We plug these solutions into equation to obtain the corresponding moment problem 


= bn [ f{T - T)e'^"'^dr 
Jo 


(47) 


for all n G N where 




(—l)^A;7r, n = 2k 
—/ife, n = 2k — 1 


(48) 


and by Proposition l2.1.21 a„ = G Again, it is easy to show that there exists constants K,6 > 0 

such that (1551) holds. From equations m and m we see that the series ^ 1/A„ converges, and that there 
exists a constant p > 0 such that |Afe+i — Afc| > p for all fc G N. This implies the existence of a biorthogonal 
sequence {6(j(T)}jgN such that there are constants Mi, M 2 > 0 such that 


||6»j|| < Mie^^^. 


Hence, as earlier, 


/(r-T) = ^ 

few 


Qj XjT 

bj 


SjiT) 


converges and provides a solution to the moment problem (|47)) . The proof of Proposition [LOT] for the case 
of Dirichlet control ([2|) , is a direct consequence of Lemma 13.2.11 

Remark 3.2.1. The numbers bj. in (H51) and dMl), are called control input coefficients and can be viewed as 
Fourier coefficients of an element b 0 /A_i for which the control problem (ED with either © or can be 
formulated as 


y = Ay + bf, y{0) = y°. 

In the Neumann case, the input element is admissible on the state space Xq =%, or equivalently that (j28l) 
and dMD hold. In the Dirichlet case, b is admissible on the state space A_i/ 2 . Both of these spaces are slightly 
suboptimal in the sense that the Carleson measure criterion due to Ho and Russell m and Weiss { 73 / can 
he used as in m to show admissibility holds in the spaces A 1/4 and X_i /4 respectively for the Neumann and 
Dirichlet control problems. 
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